Using ancient sedimentary DNA to forecast ecosystem trajectories under climate change.

Ecosystem response to climate change is complex. In order to forecast ecosystem dynamics, we need high-quality data on changes in past species abundance that can inform process-based models. Sedimentary ancient DNA (sedaDNA) has revolutionised our ability to document past ecosystems' dynamics. It provides time series of increased taxonomic resolution compared to microfossils (pollen, spores), and can often give species-level information, especially for past vascular plant and mammal abundances. Time series are much richer in information than contemporary spatial distribution information, which have been traditionally used to train models for predicting biodiversity and ecosystem responses to climate change. Here, we outline the potential contribution of sedaDNA to forecast ecosystem changes. We showcase how species-level time series may allow quantification of the effect of biotic interactions in ecosystem dynamics, and be used to estimate dispersal rates when a dense network of sites is available. By combining palaeo-time series, process-based models, and inverse modelling, we can recover the biotic and abiotic processes underlying ecosystem dynamics, which are traditionally very challenging to characterise. Dynamic models informed by sedaDNA can further be used to extrapolate beyond current dynamics and provide robust forecasts of ecosystem responses to future climate change. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

LocoGSE, a sequence-based genome size estimator for plants.

Extensive research has focused on exploring the range of genome sizes in eukaryotes, with a particular emphasis on land plants, where significant variability has been observed. Accurate estimation of genome size is essential for various research purposes, but existing sequence-based methods have limitations, particularly for low-coverage datasets. In this study, we introduce LocoGSE, a novel genome size estimator designed specifically for low-coverage datasets generated by genome skimming approaches. LocoGSE relies on mapping the reads on single copy consensus proteins without the need for a reference genome assembly. We calibrated LocoGSE using 430 low-coverage Angiosperm genome skimming datasets and compared its performance against other estimators. Our results demonstrate that LocoGSE accurately predicts monoploid genome size even at very low depth of coverage (<1X) and on highly heterozygous samples. Additionally, LocoGSE provides stable estimates across individuals with varying ploidy levels. LocoGSE fills a gap in sequence-based plant genome size estimation by offering a user-friendly and reliable tool that does not rely on high coverage or reference assemblies. We anticipate that LocoGSE will facilitate plant genome size analysis and contribute to evolutionary and ecological studies in the field. Furthermore, at the cost of an initial calibration, LocoGSE can be used in other lineages.

High resolution ancient sedimentary DNA shows that alpine plant diversity is associated with human land use and climate change.

The European Alps are highly rich in species, but their future may be threatened by ongoing changes in human land use and climate. Here, we reconstructed vegetation, temperature, human impact and livestock over the past ~12,000 years from Lake Sulsseewli, based on sedimentary ancient plant and mammal DNA, pollen, spores, chironomids, and microcharcoal. We assembled a highly-complete local DNA reference library (PhyloAlps, 3923 plant taxa), and used this to obtain an exceptionally rich sedaDNA record of 366 plant taxa. Vegetation mainly responded to climate during the early Holocene, while human activity had an additional influence on vegetation from 6 ka onwards. Land-use shifted from episodic grazing during the Neolithic and Bronze Age to agropastoralism in the Middle Ages. Associated human deforestation allowed the coexistence of plant species typically found at different elevational belts, leading to levels of plant richness that characterise the current high diversity of this region. Our findings indicate a positive association between low intensity agropastoral activities and precipitation with the maintenance of the unique subalpine and alpine plant diversity of the European Alps.

Postglacial species arrival and diversity buildup of northern ecosystems took millennia.

What drives ecosystem buildup, diversity, and stability? We assess species arrival and ecosystem changes across 16 millennia by combining regional-scale plant sedimentary ancient DNA from Fennoscandia with near-complete DNA and trait databases. We show that postglacial arrival time varies within and between plant growth forms. Further, arrival times were mainly predicted by adaptation to temperature, disturbance, and light. Major break points in ecological trait diversity were seen between 13.9 and 10.8 calibrated thousand years before the present (cal ka BP), as well as break point in functional diversity at 12.0 cal ka BP, shifting from a state of ecosystem buildup to a state where most habitat types and biotic ecosystem components were in place. Trait and functional diversity stabilized around 8 cal ka BP, after which both remained stable, although changes in climate took place and species inflow continued. Our ecosystem reconstruction indicates a millennial-scale time phase of formation to reach stable and resilient levels of diversity and functioning.

ORTHOSKIM: In silico sequence capture from genomic and transcriptomic libraries for phylogenomic and barcoding applications.

Low-coverage whole genome shotgun sequencing (or genome skimming) has emerged as a cost-effective method for acquiring genomic data in nonmodel organisms. This method provides sequence information on chloroplast genome (cpDNA), mitochondrial genome (mtDNA) and nuclear ribosomal regions (rDNA), which are over-represented within cells. However, numerous bioinformatic challenges remain to accurately and rapidly obtain such data in organisms with complex genomic structures and rearrangements, in particular for mtDNA in plants or for cpDNA in some plant families. Here we introduce the pipeline ORTHOSKIM, which performs in silico capture of targeted sequences from genomic and transcriptomic libraries without assembling whole organelle genomes. ORTHOSKIM proceeds in three steps: (i) global sequence assembly, (ii) mapping against reference sequences and (iii) target sequence extraction; importantly it also includes a range of quality control tests. Different modes are implemented to capture both coding and noncoding regions of cpDNA, mtDNA and rDNA sequences, along with predefined nuclear sequences (e.g., ultraconserved elements) or collections of single-copy orthologue genes. Moreover, aligned DNA matrices are produced for phylogenetic reconstructions, by performing multiple alignments of the captured sequences. While ORTHOSKIM is suitable for any eukaryote, a case study is presented here, using 114 genome-skimming libraries and four RNA sequencing libraries obtained for two plant families, Primulaceae and Ericaceae, the latter being a well-known problematic family for cpDNA assemblies. ORTHOSKIM recovered with high success rates cpDNA, mtDNA and rDNA sequences, well suited to accurately infer evolutionary relationships within these families. ORTHOSKIM is released under a GPL-3 licence and is available at: https://github.com/cpouchon/ORTHOSKIM.

Late Quaternary dynamics of Arctic biota from ancient environmental genomics.

During the last glacial-interglacial cycle, Arctic biotas experienced substantial climatic changes, yet the nature, extent and rate of their responses are not fully understood1-8. Here we report a large-scale environmental DNA metagenomic study of ancient plant and mammal communities, analysing 535 permafrost and lake sediment samples from across the Arctic spanning the past 50,000 years. Furthermore, we present 1,541 contemporary plant genome assemblies that were generated as reference sequences. Our study provides several insights into the long-term dynamics of the Arctic biota at the circumpolar and regional scales. Our key findings include: (1) a relatively homogeneous steppe-tundra flora dominated the Arctic during the Last Glacial Maximum, followed by regional divergence of vegetation during the Holocene epoch; (2) certain grazing animals consistently co-occurred in space and time; (3) humans appear to have been a minor factor in driving animal distributions; (4) higher effective precipitation, as well as an increase in the proportion of wetland plants, show negative effects on animal diversity; (5) the persistence of the steppe-tundra vegetation in northern Siberia enabled the late survival of several now-extinct megafauna species, including the woolly mammoth until 3.9 ± 0.2 thousand years ago (ka) and the woolly rhinoceros until 9.8 ± 0.2 ka; and (6) phylogenetic analysis of mammoth environmental DNA reveals a previously unsampled mitochondrial lineage. Our findings highlight the power of ancient environmental metagenomics analyses to advance understanding of population histories and long-term ecological dynamics.

Environmental palaeogenomic reconstruction of an Ice Age algal population.

Palaeogenomics has greatly increased our knowledge of past evolutionary and ecological change, but has been restricted to the study of species that preserve either as or within fossils. Here we show the potential of shotgun metagenomics to reveal population genomic information for a taxon that does not preserve in the body fossil record, the algae Nannochloropsis. We shotgun sequenced two lake sediment samples dated to the Last Glacial Maximum and reconstructed full chloroplast and mitochondrial genomes to explore within-lake population genomic variation. This revealed two major haplogroups for each organellar genome, which could be assigned to known varieties of N. limnetica, although we show that at least three haplotypes were present using our minimum haplotype diversity estimation method. These approaches demonstrate the utility of lake sedimentary ancient DNA (sedaDNA) for population genomic analysis, thereby opening the door to environmental palaeogenomics, which will unlock the full potential of sedaDNA.

The Treasure Vault Can be Opened: Large-Scale Genome Skimming Works Well Using Herbarium and Silica Gel Dried Material.

Genome skimming has the potential for generating large data sets for DNA barcoding and wider biodiversity genomic studies, particularly via the assembly and annotation of full chloroplast (cpDNA) and nuclear ribosomal DNA (nrDNA) sequences. We compare the success of genome skims of 2051 herbarium specimens from Norway/Polar regions with 4604 freshly collected, silica gel dried specimens mainly from the European Alps and the Carpathians. Overall, we were able to assemble the full chloroplast genome for 67% of the samples and the full nrDNA cluster for 86%. Average insert length, cover and full cpDNA and rDNA assembly were considerably higher for silica gel dried than herbarium-preserved material. However, complete plastid genomes were still assembled for 54% of herbarium samples compared to 70% of silica dried samples. Moreover, there was comparable recovery of coding genes from both tissue sources (121 for silica gel dried and 118 for herbarium material) and only minor differences in assembly success of standard barcodes between silica dried (89% ITS2, 96% matK and rbcL) and herbarium material (87% ITS2, 98% matK and rbcL). The success rate was > 90% for all three markers in 1034 of 1036 genera in 160 families, and only Boraginaceae worked poorly, with 7 genera failing. Our study shows that large-scale genome skims are feasible and work well across most of the land plant families and genera we tested, independently of material type. It is therefore an efficient method for increasing the availability of plant biodiversity genomic data to support a multitude of downstream applications.

Plant DNA metabarcoding of lake sediments: How does it represent the contemporary vegetation.

Metabarcoding of lake sediments have been shown to reveal current and past biodiversity, but little is known about the degree to which taxa growing in the vegetation are represented in environmental DNA (eDNA) records. We analysed composition of lake and catchment vegetation and vascular plant eDNA at 11 lakes in northern Norway. Out of 489 records of taxa growing within 2 m from the lake shore, 17-49% (mean 31%) of the identifiable taxa recorded were detected with eDNA. Of the 217 eDNA records of 47 plant taxa in the 11 lakes, 73% and 12% matched taxa recorded in vegetation surveys within 2 m and up to about 50 m away from the lakeshore, respectively, whereas 16% were not recorded in the vegetation surveys of the same lake. The latter include taxa likely overlooked in the vegetation surveys or growing outside the survey area. The percentages detected were 61, 47, 25, and 15 for dominant, common, scattered, and rare taxa, respectively. Similar numbers for aquatic plants were 88, 88, 33 and 62%, respectively. Detection rate and taxonomic resolution varied among plant families and functional groups with good detection of e.g. Ericaceae, Roseaceae, deciduous trees, ferns, club mosses and aquatics. The representation of terrestrial taxa in eDNA depends on both their distance from the sampling site and their abundance and is sufficient for recording vegetation types. For aquatic vegetation, eDNA may be comparable with, or even superior to, in-lake vegetation surveys and may therefore be used as an tool for biomonitoring. For reconstruction of terrestrial vegetation, technical improvements and more intensive sampling is needed to detect a higher proportion of rare taxa although DNA of some taxa may never reach the lake sediments due to taphonomical constrains. Nevertheless, eDNA performs similar to conventional methods of pollen and macrofossil analyses and may therefore be an important tool for reconstruction of past vegetation.

Ancient plant DNA in lake sediments.

Contents 924 I. 925 II. 925 III. 927 IV. 929 V. 930 VI. 930 VII. 931 VIII. 933 IX. 935 X. 936 XI. 938 938 References 938 SUMMARY: Recent advances in sequencing technologies now permit the analyses of plant DNA from fossil samples (ancient plant DNA, plant aDNA), and thus enable the molecular reconstruction of palaeofloras. Hitherto, ancient frozen soils have proved excellent in preserving DNA molecules, and have thus been the most commonly used source of plant aDNA. However, DNA from soil mainly represents taxa growing a few metres from the sampling point. Lakes have larger catchment areas and recent studies have suggested that plant aDNA from lake sediments is a more powerful tool for palaeofloristic reconstruction. Furthermore, lakes can be found globally in nearly all environments, and are therefore not limited to perennially frozen areas. Here, we review the latest approaches and methods for the study of plant aDNA from lake sediments and discuss the progress made up to the present. We argue that aDNA analyses add new and additional perspectives for the study of ancient plant populations and, in time, will provide higher taxonomic resolution and more precise estimation of abundance. Despite this, key questions and challenges remain for such plant aDNA studies. Finally, we provide guidelines on technical issues, including lake selection, and we suggest directions for future research on plant aDNA studies in lake sediments.

Lake sedimentary DNA accurately records 20th Century introductions of exotic conifers in Scotland.

Sedimentary DNA (sedDNA) has recently emerged as a new proxy for reconstructing past vegetation, but its taphonomy, source area and representation biases need better assessment. We investigated how sedDNA in recent sediments of two small Scottish lakes reflects a major vegetation change, using well-documented 20th Century plantations of exotic conifers as an experimental system. We used next-generation sequencing to barcode sedDNA retrieved from subrecent lake sediments. For comparison, pollen was analysed from the same samples. The sedDNA record contains 73 taxa (mainly genus or species), all but one of which are present in the study area. Pollen and sedDNA shared 35% of taxa, which partly reflects a difference in source area. More aquatic taxa were recorded in sedDNA, whereas taxa assumed to be of regional rather than local origin were recorded only as pollen. The chronology of the sediments and planting records are well aligned, and sedDNA of exotic conifers appears in high quantities with the establishment of plantations around the lakes. SedDNA recorded other changes in local vegetation that accompanied afforestation. There were no signs of DNA leaching in the sediments or DNA originating from pollen.

The role of sea ice for vascular plant dispersal in the Arctic.

Sea ice has been suggested to be an important factor for dispersal of vascular plants in the Arctic. To assess its role for postglacial colonization in the North Atlantic region, we compiled data on the first Late Glacial to Holocene occurrence of vascular plant species in East Greenland, Iceland, the Faroe Islands and Svalbard. For each record, we reconstructed likely past dispersal events using data on species distributions and genetics. We compared these data to sea-ice reconstructions to evaluate the potential role of sea ice in these past colonization events and finally evaluated these results using a compilation of driftwood records as an independent source of evidence that sea ice can disperse biological material. Our results show that sea ice was, in general, more prevalent along the most likely dispersal routes at times of assumed first colonization than along other possible routes. Also, driftwood is frequently dispersed in regions that have sea ice today. Thus, sea ice may act as an important dispersal agent. Melting sea ice may hamper future dispersal of Arctic plants and thereby cause more genetic differentiation. It may also limit the northwards expansion of competing boreal species, and hence favour the persistence of Arctic species.

Characterization of 14 microsatellite markers for Silene acaulis (Caryophyllaceae)(1).

PREMISE OF THE STUDY: Fifty candidate microsatellite markers, generated using 454 shotgun sequencing, were tested for the widespread arctic/alpine herb Silene acaulis (Caryophyllaceae). METHODS AND RESULTS: Fourteen out of 50 markers resulted in polymorphic products with profiles that enabled interpretation. The numbers of alleles per locus ranged from two to six, and the expected heterozygosity per locus ranged from 0.06 to 0.68. Analysis of F0 and F1 samples proved that one allele was always inherited maternally. Four multiplex mixes have been developed. CONCLUSIONS: Microsatellite markers for this species will be a valuable tool to study detailed small-scale genetic patterns in an arctic/alpine herb and to relate them to demographic parameters.

Comparative analyses of plastid and AFLP data suggest different colonization history and asymmetric hybridization between Betula pubescens and B. nana.

Birches (Betula spp.) hybridize readily, confounding genetic signatures of refugial isolation and postglacial migration. We aimed to distinguish hybridization from range-shift processes in the two widespread and cold-adapted species B. nana and B. pubescens, previously shown to share a similarly east-west-structured variation in plastid DNA (pDNA). We sampled the two species throughout their ranges and included reference samples of five other Betula species and putative hybrids. We analysed 901 individual plants using mainly nuclear high-resolution markers (amplified fragment length polymorphisms; AFLPs); a subset of 64 plants was also sequenced for two pDNA regions. Whereas the pDNA variation as expected was largely shared between B. nana and B. pubescens, the two species were distinctly differentiated at AFLP loci. In B. nana, both the AFLP and pDNA results corroborated the former pDNA-based hypothesis that it expanded from at least two major refugia in Eurasia, one south of and one east of the North European ice sheets. In contrast, B. pubescens showed a striking lack of geographic structuring of its AFLP variation. We identified a weak but significant increase in nuclear (AFLP) gene flow from B. nana into B. pubescens with increasing latitude, suggesting hybridization has been most frequent at the postglacial expansion front of B. pubescens and that hybrids mainly backcrossed to B. pubescens. Incongruence between pDNA and AFLP variation in B. pubescens can be explained by efficient expansion from a single large refugium combined with leading-edge hybridization and plastid capture from B. nana during colonization of new territory already occupied by this more cold-tolerant species.

Long-distance plant dispersal to North Atlantic islands: colonization routes and founder effect.

Long-distance dispersal (LDD) processes influence the founder effect on islands. We use genetic data for 25 Atlantic species and similarities among regional floras to analyse colonization, and test whether the genetic founder effect on five islands is associated with dispersal distance, island size and species traits. Most species colonized postglacially via multiple dispersal events from several source regions situated 280 to >3000 km away, and often not from the closest ones. A strong founder effect was observed for insect-pollinated mixed maters, and it increased with dispersal distance and decreased with island size in accordance with the theory of island biogeography. Only a minor founder effect was observed for wind-pollinated outcrossing species. Colonization patterns were largely congruent, indicating that despite the importance of stochasticity, LDD is mainly determined by common factors, probably dispersal vectors. Our findings caution against a priori assuming a single, close source region in biogeographic analyses.

Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

Genetic roadmap of the Arctic: plant dispersal highways, traffic barriers and capitals of diversity.

We provide the first comparative multispecies analysis of spatial genetic structure and diversity in the circumpolar Arctic using a common strategy for sampling and genetic analyses. We aimed to identify and explain potential general patterns of genetic discontinuity/connectivity and diversity, and to compare our findings with previously published hypotheses. We collected and analyzed 7707 samples of 17 widespread arctic-alpine plant species for amplified fragment length polymorphisms (AFLPs). Genetic structure, diversity and distinctiveness were analyzed for each species, and extrapolated to cover the geographic range of each species. The resulting maps were overlaid to produce metamaps. The Arctic and Atlantic Oceans, the Greenlandic ice cap, the Urals, and lowland areas between southern mountain ranges and the Arctic were the strongest barriers against gene flow. Diversity was highest in Beringia and gradually decreased into formerly glaciated areas. The highest degrees of distinctiveness were observed in Siberia. We conclude that large-scale general patterns exist in the Arctic, shaped by the Pleistocene glaciations combined with long-standing physical barriers against gene flow. Beringia served as both refugium and source for interglacial (re)colonization, whereas areas further west in Siberia served as refugia, but less as sources for (re)colonization.